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ABSTRACT 

We present a near-infrared quasar composite spectrum spanning the wave- 
length range 0.58 - 3.5 /im. The spectrum has been constructed from obser- 
vations of 27 quasars obtained at the NASA IRTF telescope and satisfying the 
criteria Kg < 14.5 and < -23; the redshift range is 0.118 < z < 0.418. The 
signal-to-noise is moderate, reaching a maximum of 150 between 1.6 and 1.9 /im. 
While a power-law fit to the continuum of the composite spectrum requires two 
breaks, a single power-law slope of a = —0.92 plus a 1260 K blackbody provides 
an excellent description of the spectrum from Ha to 3.5 fim, strongly suggest- 
ing the presence of significant quantities of hot dust in this blue-selected quasar 
sample. We measure intensities and line widths for ten lines, finding that the 
Paschen line ratios rule out Case B recombination. We compute i^-corrections 
for the J, H, K, and Spitzer 3.6 /im bands which will be useful in analyzing 
observations of quasars up to 2; = 10. 

Subject headings: quasars: emission lines — quasars: general 



^Columbia Univeristy Department of Astronomy, 550 West 120th Street, New York, NY 10027 
^Space Telescope Science Institute, 3700 San Martin Dr., Baltimore, MD 21218 



- 2 - 



1. Introduction 

Spectral templates constructed from large uniform source samples have been invaluable 
for understanding the class of objects in question. Composite spectra enhance common 
features, while unusual characteristics are averaged out (see discussions in Scott et al. 2004; 
Telfer et al. 2002; Vanden Berk et al. 2001; Brotherton et al. 2001; Zheng et al. 1997; Francis 
et al. 1991; Boyle 1990; Cristiani & Vic 1990). For quasars in the ultraviolet to optical 
regimes, comparison of individual spectra with such composites helps to distinguish normal 
features from those which arise from unusual physical conditions. The average shape of 
the rest-frame UV continuum has provided insight into the source of photoionization of the 
quasar's emission line gas, while the emission lines provide information on temperatures and 
densities and have helped discriminate among AGN models. 

Composite spectra can also be useful in quasar searches, and in comparing new sub- 
classes of sources with existing populations. For example, we (Clikman et al. 2004) have 
recently sought to assess the incompleteness of the blue-selected quasar population by search- 
ing for highly obscured and/or very red quasars among a candidate list derived by matching 
the FIRST (radio) and 2MASS (near-IR) surveys. We found seventeen highly obscured 
quasars (E(B-V)> 0.5) with red spectral energy distributions by obtaining near-IR spectra 
spanning the wavelength range 0.9-2.5//m. Since existing quasar composite spectra extend 
only to ~ 8500 A, they are of little use in interpreting the line and continuum emission in 
these obscured objects for z < 1. This motivated us to construct a near-infrared composite 
spectrum for a sample of normal quasars with blue colors and little or no obscuration in 
order to be able to assess how similar our red-selected sample is to the dominant population. 
In addition, our IR composite will prove useful when observing high-redshift objects with 
the InfraRed Spectrograph (IRS) onboard Spitzer whose wavelength coverage extends from 
5 to 35/xm, allowing comparisons for objects out to 2; ~ 4 — 5. 

In this work we present the first near-infrared quasar composite spectrum covering the 
0.58-3.5;Lim rest wavelength range. Our quasars were chosen to be luminous, low-redshift, 
UVX-selected quasars from the Sloan Digital Sky Survey Quasar Catalog First Data Release 
(SDSS-DRl) catalog. We discuss the details of our selection criteria and the observations 
in §2. We describe the method for constructing the composite and present the results in §3. 
In §4 we analyze the continuum and line features in the spectrum followed by a comparison 
with previous infrared quasar SEDs from the literature (§5). In §6 we present a combined 
optical-to-near-infrared composite, utilizing the SDSS optical spectrum available for each 
source, and use the resultant spectrum to compute /^-corrections for the 2MASS J, H , and 
K bands, as well as for the Spitzer 3.6 /xm IRAC band. We summarize this work in §7 and 
discuss future efforts to extend and improve our knowledge of the intrinsic quasar spectrum. 
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We use standard cosmological parameters throughout the paper: Hq — 70 km s ^ Mpc ^, 
Qm = 0.30, Qa = 0.70. 

2. Object Selection and Observations 

The quasars selected for inclusion in this composite spectrum were chosen from the 
SDSS-DRl quasar catalog (Schneider et al. 2003) which contains 16,713 objects. Of these, 
2260 had detections in the Two Micron All Sky Survey (2MASS, Strutskie et al. 1997) Point 
Source Catalog (PSC). In order to optimize the use of telescope time by obtaining many 
spectra with good signal-to-noise ratios, we selected objects with Ks < 14.5. We restrict 
our objects to relatively low redshifts {z < 0.5) to focus on the quasars' rest-frame near- 
infrared emission; higher redshift objects will contain optical lines which have already been 
explored in optical quasar surveys. Since we would like to avoid a significant contribution 
from the quasar's host galaxy in the K band, we require that our objects have absolute 
i-band magnitudes of Mj < —23. With this choice, an L* host will contribute at most 40% 
of the total i-band flux, assuming M* — —22.03 (Blanton et al. 2001); with our projected 
slit widths (see below), the galaxy fraction is < 8%. There are 87 objects in the DRl Quasar 
catalog which meet these criteria. 

We obtained near-infrared spectra for 27 of the quasars obeying the above criteria. The 
observations are summarized in Table 1. All spectra but one were obtained at the NASA 
Infrared Telescope Facility (IRTF), a 3.0 meter telescope optimized for infrared observations 
on the summit of Mauna Kea, Hawai'i. We used Spex (Rayner et al. 2003) in short and 
long cross-dispersed modes (SXD, 0.8-2.5 fim and LXD, 1.9-4.1 ^m, respectively), obtaining 
24 SXD and 13 LXD spectra with A/(AA) ~ 1200 and 1500, respectively The data in 
each mode were reduced using Spextool (Gushing, Vacca, & Rayner 2004) and corrected 
for telluric absorption and instrument response with the Telluric software package (Vacca, 
Gushing, & Rayner 2003); a standard AOV star spectrum was obtained immediately after 
each object spectrum using calibrators that lay within both a difference in airmass of 0.1 
and a separation from the target of < 10°. The SXD and LXD spectra were merged using 
the MERGEXD task within the Spextool package which scales the weaker spectrum to the 
brighter spectrum using a selected region in the K band where the two modes overlap. 
The fully reduced spectra are available in the online version of the article. The IRTF data 
were obtained between 2003 September 23 and 2004 September 14 over three observing 
runs. Seeing ranged from ~ 0'.'6 to ~ 1'.'25 FWHM; we mostly used the 0'.'5 and 0'.'8 slits as 
appropriate (see columns (6) and (9) of Table 1). For all but a few of the SXD spectra, we 
used the 0'/5 slit which projects 1.5 — 2.8 kpc on the sky in our redshift range 0.18 — 0.42. This 
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is at least an order of magnitude larger than the size of the typical narrow-line region of a 
quasar (~ 100 pc) and wc therefore are inevitably including in a large amount of host galaxy 
light. Assuming a de Vaucouleurs (r^/'^-law) surface-density profile with a half-light radius 
of 11.5 kpc (Dunlop et al. 2003) and depending on the quasar's redshift, this corresponds 
to 13 — 20% of the total host galaxy light. We discuss in greater detail the impact of host 
galaxy light on our composite in Sections 4.1 and 6. 

One spectrum was obtained at the MDM Observatory on Kitt Peak on 2004 November 
24 with the TIFKAM infrared camera in spectroscopic mode using a 0'.'6 slit . Conditions 
were photometric and the seeing was 0'.'9 FWHM. This spectrum spans 0.95-2.45/im at a 
lower resolving power of 600. The spectrum was extracted and reduced using standard 
spectral reduction tasks in IRAF. The telluric absorption and instrument response were 
corrected using a general version of the Telluric software written to support non-SpeX data. 

The 27 observed sources and some of their known properties are listed in Table 2. 
Columns (1) and (2) hst the quasars' J2000.0 coordinates from the SDSS-DRl catalog. 
Column (3) hsts the redshift from the SDSS quasar catalog, followed by the SDSS u,g,r,i, 
and z magnitudes in columns (4) - (8). These are on the AB magnitude system defined 
by ASmag = -2.5 logio(/i./3631Jy). The 2MASS J, H, and K magnitudes are hsted in 
columns (9), (10), and (11), respectively; these arc Vega-based, not AB, magnitudes. About 
one-third of the objects wc observed have detections in the FIRST survey (Becker. White. 
& Helfand 1995) implying 5*20^1 > 1-0 mJy; columns (12) and (13) list the 20 cm peak flux 
density from the FIRST survey (wc report the NVSS value (Condon ct al. 1998) for one 
object which falls outside the FIRST survey area) and the radio-loudness parameter (the 
ratio of radio to optical flux density Rr = fraxHo/ f optical) with respect to the r band. Column 
(14) reports the absolute i-band magnitude, Mj. 



3. Constructing the Composite 

When the extraction procedure is applied to the different orders of the cross-dispersed 
IRTF spectra, it yields spectra with non-linear dispersions. In the /-band order spanning 
the wavelengths 0.83 - 0.93/xm, the dispersion is about 2.3A/pixel and is roughly constant 
until the z-band (0.93 - 1.06//m) where the dispersion is about 2.7A/pixel. The J, H, and 
K orders of the SXD mode have wavelength increments of about 3.2, 4.0, and 5.lA/pixel, 
respectively. The K-band orders (1.9 - 2.55/im and 2.55 - 3.14yum) in the LXD mode have 
larger dispersions of about 5.4 and 5.8A while the two L band orders (2.98 - 3.52//m and 
3.52 - 4.2/im) have wavelength increments of about 6.8 and 8.1 A/pixel. These values are 
relatively constant within the spectral orders, varying by ±0.2A/pixel per order. 
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To create a mean spectrum it is necessary to shift the objects to their rest frame using 
the redshifts from the SDSS quasar catalog and to rebin everything to a common dispersion. 
This causes the rest frame orders to overlap, presenting a challenge for choosing a common 
dispersion that preserves the most information without oversampling in regions of lower dis- 
persion. Although the dispersion changes by more than a factor of four over the wavelength 
range of our spectra, the resolving power, TZ = A/(AA), only changes from 3700 to 5000. 
Therefore, we chose to rebin our spectra by picking a dispersion that is proportional to the 
wavelength so that 7?. is a constant. With this choice, shifting spectra at different redshifts 
will not change TZ. We chose to fix TZ to its median value in the observed K band since all 
our spectra cover this region. We measure (A = 21930A)/(AA = 5.lA) = 4300 and build a 
wavelength grid using this dispersion in log A increments starting at the shortest rest-frame 
wavelength in our data, Ao = 5710A. The recursion relation for the wavelength is 



" (1 - 1/4300)'^ ■ ^ ^ 

The longest rest-frame wavelength in our spectra is 37000A, which requires roughly 
n — 8030 elements in our wavelength grid. In this wavelength grid, the dispersion ranges 
smoothly from about 2A to lOA. We linearly interpolate the flux and error arrays onto this 
wavelength grid. 

We built our composite spectrum following the method outhned in Francis et al. (1991). 

We ordered the quasars by increasing redshift and scaled the lowest redshift spectrum by 
its median value, excluding regions of atmospheric absorption. After removing atmospheric 
absorption regions and shifting to the rest frame, wc then scaled the next spectrum to the 
first object using the median fiux in their overlapping wavelength regions. By using the 
median fiux for scaling each subsequent spectrum, we avoid the influence of strong emission 
lines. We weight each spectrum by the inverse square of its error spectrum and sum these 
products. We scale subsequent spectra to the median flux in the region overlapping the 
partial composite from the previous cycle. The flnal weighted sum is divided by the total of 
the weights to yield a weighted average composite spectrum: 



<h>- ^n^y^j^ . {2) 

Finally, we smooth the resulting spectrum with a boxcar, averaging over 5 pixels. Figure 
1 displays the flnal composite with the expected prominent emission lines highlighted by 
dotted lines. 



The upper panel of Figure 2 displays the corresponding combined error spectrum (on 
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an expanded y-axis) obtained by propagating the individual error spectra according to the 
equation, 



and smoothing in the same manner as for the composite spectrum. We compute the signal- 
to-noise ratio per pixel by dividing the composite by the combined error spectrum. The 
lower panel of Figure 2 shows the signal-to-noise as a function of wavelength. We achieve 
a signal-to-noise ratio of ~ 150 in the 1.6-1.9 /im wavelength region, which, in most of our 
objects, corresponds to the observed K-hand. This is where the most photons were collected 
as a consequence of the overlap of the two SpeX modes in this region. Spectra covering the 
observed 0.8 - 4.1/im range have had the overlapping i^-band region observed twice; the 
duplicate observations were combined to create the final object spectrum. 

Assuming the quasar continuum is a power-law, fx oc X~(°''^+'^\ a geometric mean spec- 
trum enables a power-law fit to the continuum that is representative of the arithmetic mean 
of the indices of the individual spectra, < fx >gmOC X-^<°'''>+'^\ We have also generated 
a geometric mean spectrum, < fx >gm, which does not apply a weighting from the error 
spectra, according to the following formalism: 



where n is the number of spectra contributing to the ith wavelength unit. We scale the 
spectra using the same scaling as we did for the weighted arithmetic mean, excluding regions 
of atmospheric absorption in each spectrum. Figure 3 displays the geometric mean spectrum 
marked with the same emission lines as in Figure 1. Figure 4 plots the number of spectra 
contributing at each wavelength to the geometric mean spectrum. Table 3 lists the arithmetic 
mean spectrum, its error spectrum, and the geometric mean spectrum as a function of 
wavelength. 

The order in which spectra are added affects the final composite spectrum. This was 
noted by Brotherton et al. (2001) who found a change in the spectral index of their composite 
from a = —0.46 {F^, oc z/") when combining the spectra starting from the lowest rcdshift 
to a = —0.49 when starting with the highest rcdshift quasars. The first few spectra in 
the composite-building algorithm can strongly influence the shape of the final composite 
if they contain wavelength regions that exist in only a few spectra. Telfer et al. (2002) 
order their spectra beginning with objects that contribute to the best-sampled part of their 
wavelength coverage. We constructed two additional composite quasar spectra to study the 




(3) 



n 



</A>gm= (11 /aO'/", 
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dependence of the spectral shape of our composite on the order in which we combine the 
spectra. Following Brotherton et al. (2001), we made one spectrum starting with the highest 
redshift quasar and adding sequently lower redshift objects to the composite. To make a 
composite using the method of Telfer et al. (2002), we detemined (as shown in Figure 2) 
that quasars within the redshift range 0.188 < z < 0.316 contribute to the 1.6 — 1.9 /xm 
wavelength region, which has the highest signal-to-noise in the composite. We started with 
SDSSJ015910.0+010514.5, a z = 0.2174 quasar that has both SXD and LXD coverage and 
added to it spectra that have SXD and LXD coverage in same wavelength range. We then 
added spectra in that wavelength range having only an SXD or LXD spectrum. Finally, we 
added sources with redshifts outside this high-SNR range, alternating between lower and 
higher redshift objects. We discuss the sensitivity of the resultant spectra to these different 
ordering methods in §4. The analysis of our near-infrared quasar composite described below 
pertains to the composite made from quasars added in increasing redshift order. 

Wc also constructed a uniformly weighted composite spectrum in order to verify that our 
weighting scheme does not skew our composite toward the highest signal-to-noise spectra. 
In this case, instead of using equation (2) in our composite-building algorithm, we use 
< fx >= 'Y^^i f\i/n^ excluding the same regions of atmospheric absorption as in the error- 
weighted composite. The uniformly weighted spectrum has a similar continuum shape and 
line profiles as the error- weighted composite. We verified this by examining the ratio of, and 
difference between, the two spectra. No obvious feature emerged and there were no residuals 
at the location of emission lines. 

The signal-to-noise ratios in the observed i^-band for the individual spectra that make 
up our composite range from ~ 1.3 — 60; the median signal-to-noise ratio is 12.4, the mean 
is 16, and fourteen of the 27 quasars in our sample have signal-to-noise ratios between 10 
and 30. Comparison of our uniformly weighted composite and the highest signal-to-noise 
spectrum shows strong residuals both in the continuum and at the position of emission lines. 
We conclude that our error-weighted composite is not dominated by the highest signal-to- 
noise spectra in our sample. 



4. Analysis of the Composite Spectrum 

4.1. Is Our Composite Biased? 

To evaluate any bias our selection criteria may introduce with respect to the parent 
quasar sample, wc constructed an optical composite spectrum using the technique described 
above from the 27 SDSS optical spectra available for our sample. We rebinned the spectra 
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using the same resolving power, TZ — 4300, and applied the same scaling technique after 
sorting the spectra by increasing redshift. We also produced optical composites using the 
two additional sorting methods: we sorted the quasars by decreasing redshift and sorted 
according to the objects whose redshifts contribute to the highest signal-to-noise part of 
the final composite in the same manner described in §3. The SDSS database provides an 
error spectrum for each object which we use to weight the data in the same manner as for 
the near-infrared composite. The top panel of Figure 5 plots the resultant optical composite 
(black) for our 27 objects along with a composite that we generated from the complete sample 
of 87 SDSS-DRl quasars that obeyed our selection criteria (blue). We compare these with 
the SDSS quasar composite spectrum created from the parent DRl sample by Vanden Berk 
et al. (2001) (red). Our selection criteria do introduce a small bias which leads to a quasar 
with a bluer continuum than the total SDSS sample composite. 

Vanden Berk et al. (2001) fit two power laws to the quasar continuum, fitting one spectral 
index for wavelengths shorter than 4000 A and one for wavelengths longer than 4000 A. To 
compute the spectral index for wavelengths longer than 4000 A, Vanden Berk et al. (2001) 
fit a power law to the line-free wavelength ranges 6006-6035 A and 7160-7180 A. They find 
a spectral index of = —2.45 using their geometric mean spectrum and a^, = —1.58 using 
their median-combined spectrum. We fit a power law to the same wavelength regions in 
our optical composite and measure a spectral index of ai, — —1.03 ± 0.1. In Section 4.2 
we analyze the continuum properties of our near-infrared composite spectrum which also 
measures the spectral index around Ha. These measurements are shown in Table 4. The 
wavelength regions used by Vanden Berk et al. (2001) to fit the blue end of their composite 
(1350-1365 and 4200-4230 A) are not available using our composite made of quasars with 
z < 0.5. However it is clear from Figure 5 that the SDSS composite has a softer UV spectrum 
than that produced from our subsample. 

This discrepancy is likely a consequence of the varying properties of quasars that con- 
tribute to the blue and red ends of the SDSS composite spectrum^. As the composite moves 
to longer rest-wavelengths, more lower-^ objects with lower average luminosities contribute. 
The median absolute i magnitude, Mj, for our sample is —23.7 and our maximum redshift 
is 2; = 0.418. In contrast, the objects contributing to the blue end of the SDSS composite 
spectrum can have redshifts as high as 2; ~ 1.7, and the median Mj for objects with z < 1.7 is 
—25.0. At the red end of the SDSS spectrum, however, only objects with z < 0.4 contribute; 
their median Mj is —22.5, fainter than our requirement of Mi < —23.0. 



^We find no distinction in the red to near-infrared {i — Ks) colors between our sample and the SDSS-DRl 

objects with the same redshift limit as our sample (z < 0.5). The median i — colors for our sample and 
the SDSS-DRl quasars with z < 0.5 are 3.1 and 3.2, respectively. 
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The difference spectrum is illustrated in the lower panel of Figure 5. For comparison, 
we overplot the elliptical galaxy template spectrum from Mannucci et al. (2001) (green) to 
determine if any galaxy features are present in the residuals. The result strongly suggests 
that most if not all of the excess red emission in the SDSS quasar composite spectrum is 
due to host galaxy light. We identify with dotted lines obvious continuum and absorption 
line features which appear in our residual spectrum as well as the galaxy template. We 
shade the prominent TiO absorption band which appears in both spectra between 7050A 
and 7275A. We conclude that the inclusion of low-redshift, low-luminosity AGN in the SDSS 
composite on the red end of the spectrum adds host galaxy light which enhances the red 
continuum. The sources in our sample were selected to be luminous enough that host galaxy 
light would not be a significant factor in the continuum. In addition, our slit width (0'.'5 
in the SXD mode) takes in only half as much galaxy light as the SDSS spectroscopic fibers 
(3") assuming an elliptical host with a de Vaucouleurs (r^/'^-law) surface-density profile and 
a half-light radius of 11.5 kpc. 

Furthermore, in a comparison between SDSS and Palomar-Green quasars from the 
Bright Quasar Survey (BQS Schmidt & Green 1983) Jester et al. (2005) found that quasar 
samples with brighter apparent magnitudes (in B, g and i bands) are biased to the blue 
compared with a fainter sample (see their §4.1.3, and Figure 10). Our sample also selects 
bright objects because of our Kg < 14.5 requirement. In the SDSS-DRl catalog, the median 
g and i magnitudes for quasars with z < 0.5 and Mj < —23.0 but no magnitude restriction 
are 18.33 and 18.98, respectively. By imposing the Kg < 14.5 magnitude hmit, the median 
g and i magnitude in our sample are 17.27 and 16.96, respectively. Therefore the bluer 
continuum of our composite compared with the SDSS composite is also due to the color bias 
introduced by our bright magnitude criterion. 



4.2. Continuum 

A single power law does not fit the continuum from 0.58-3.5//m well. This is evident 
in Figure 6 where we display the composite in A vs. F^, space on logarithmic axes. We fit a 
broken power-law (with two breaks) to the spectrum using the IDL fitting routine AMOEBA, 
a robust fitting algorithm that performs a minimization for a given model and data while 
ignoring outliers; the algorithm does not return errors on the fitted parameters. However, 
AMOEBA does require a fitting tolerance parameter, which we specify to be 1.0 x 10~^. 
The broken power law model has six parameters: the two break points in the power law, the 
amplitude of the function at the first break and the three power law indices, one for each 
segment. We exclude Balmer, Paschen and Brackett lines as well as major ionic fines (see 
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§4.3) from the fit. Tfie wavelength windows used in our fit are marked with a thick hne at 
the bottom of Figure 6. 

Table 4 compares the spectral indices derived for different regions of our spectra. For 
the composite generated by sorting the spectra from lowest to highest redshifts (row 6), we 
find that the optical part of the spectrum containing Ha is best fit by a power law starting 
at 5700 A and breaking at 10850 A with a spectral index of a = —0.78 (where F,^ oc 
This spectral index is considerably harder than that found by (Vanden Berk et al. 2001) in 
the 6005 — 7180 A region (q;,^sdss = ~2.45), although it is consistent with the 6005 — 7180 
A slope of —1.03 derived for our optical composite (row 3). The 1.085 — 2.230 fim. region is 
best fit with an index of a = —1.81, while the 2.23-3.5 /xm region is fit by a harder spectral 
index of a = —1.03. 

It is clear from Table 4 that spectral indices measured in the UV-optical part of com- 
posite quasar spectra agree (see also Table 5 in Vanden Berk et al. 2001). The order in 
which spectra are combined does not have a strong affect on the spectral slope in the 1-2 
/im and 2-3.5 /im regions. The largest discrepancies occur in the red-optical region around 
Ha. Here, the order of combining the spectra introduces a variation of up to Aa^ ~ 0.5. 

We also note that the second break in the spectrum occurs in the noisiest part of the 
composite where the LXD order begins to contribute. This parameter is also sensitive to 
the order in which the spectra are combined, shifting the break by up to AA ~ 0.3 iim. The 
two slopes are very different and a break must occur between the two regions; however, the 
location of the break may shift when more spectra can be added to the composite to improve 
its signal-to-noise ratio. 

Excess emission between ~ 2 and 10 /im has been widely noted in previous work and has 
been attributed to hot dust at ~ 1500 K in some studies of AGN spectral energy distributions 
(SEDs) (Rieke 1978; Edelson & Malkan 1986; Barvainis 1987; Elvis ct al. 1994). We fit a 
power-law plus a blackbody to our composite spectrum with AMOEBA. We define the 
following function, 

F, = CiZ/" + C25,(Trf„,,), (5) 

where By{Td) is the Planck function, allowing AMOEBA to fit Ci, a, C2, and T^ust- We 
initialize the scaling factors to unity and the spectral slope, a, to —0.78, the value for 
the optical portion of the spectrum obtained from our broken power-law fits, and the dust 
temperature, Tdusti to 1500 K as found by Barvainis (1987). The fit converged on a spectral 
slope of —0.92 and a temperature of 1260 K. We list the result of this fit in Table 5 along 
with fits to composites generated by ordering the spectra in the different schemes described 
in §3. The order in which the spectra are combined in generating the composites does not 
effect the fitted temperature or slope by more than 10°K or Act,^ ~ 0.02. Figure 7 shows 
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the results of our fit to tfie composite ordered by increasing redshift. This model reproduces 
the continuum shape at least as well as the broken power-law fits (fits to both models yield 
~ 3) and uses two fewer parameters and a stronger physical motivation. This strongly 
suggests a hot dust component exists in the luminous blue quasars observed in our sample. 

The detection of sihcate emission at 10 and 18 //m in 5pzteer spectra of PG quasars also 
suggests dust, though at lower temperatures, ranging from 140 K to 220 K (Hao et al. 2005). 
This is not inconsistent with our result, since different dust temperatures are expected at 
different radial distances from the AGN. Barvainis (1987) fit a more complex model including 
a hot dust component over a range of radii to the SEDs of 3C273 and IRAS 13349+2438 
and found temperatures ranging from 100 K to 1500 K. The wavelength coverage of our 
composite spectrum allows us to probe the hottest grains in this range. 

4.3. Emission lines 

We manually identified emission lines in the composite spectrum and measured their in- 
tensities by fitting a Voigt profile to each isolated line. The Voigt fitting algorithm computes 
a linear continuum (2 parameters) as well as a total intensity, central peak position, doppler 
width and damping width (4 parameters). We allowed the peak position in the fit to vary. 
We integrated the continuum component of the fit and subtracted it from the integrated in- 
tensity over the wavelength region containing at least 5% of the peak fitted flux to compute 
the intensity of the line above the continuum. We also computed the equivalent width of 
each line using the fitted profile. In the case of the blended He 1 -|- Pa7 pair, we fitted a 
double Voigt profile which would require ten parameters (2 continuum parameters and two 
sets of four parameters for the individual Voigt profiles). Owing to the strong blending of 
the lines, we fixed their central peak wavelengths to the known laboratory values of 10830A 
and 1094lA, respectively. We then computed the intensity of each line with a single Voigt 
profile using the six fitted parameters pertaining to each line, also over the region containing 
> 5% of the fitted peak flux. 

The low signal-to-noise of some of the weaker lines in our spectrum made the Voigt 
fitting algorithm a better choice than a multiple-component gaussian fit. Single gaussian 
profiles do not model well the broad wings of the emsission lines and generate a poor fit. 
Additional fitting components account for these wings. However, in poor signal-to-noise lines 
(such as [S III] A 9069 A and Fe II A 9202 A) the additional gaussian components will fit 
the noise and the fitted lineshape will appear distorted. The fixed shape of the Voigt profile, 
with its four line parameters provide more realistic fits to our lines, given their range in 
signal-to-noise. 
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Figure 8 displays the emission lines over their fitted regions. The line profile is plotted 

as a dotted line over the region used to integrate the intensity and the continuum is plotted 
in a dash-dot line. Table 6 lists the line parameters extracted from the fits. Column (1) 
identifies the feature. Column (2) lists the laboratory wavelength for the feature followed by 
the fitted peak position of the line in column (3) (except in the case of the blended He I + Paj 
pair where we fixed these values). Columns (4) and (5) list the lower and upper wavelength 
limits for integrating the line intensity. Column (6) lists the line intensity relative to H-ct 
followed by the equivalent width of the hne in column (7). Column (8) lists the full- width 
at half-maximum of the line. 

The relative fluxes of Paa and Pa/3 give a line ratio of Paa/Pa/? = 0.64 ± 0.01. This 
value is less than half that predicted by Case B recombination which gives a value of 1.8 for a 
typical broad hne cloud with Tg = 10^ K and Ne — 10^ cm~^. The expected value is constant 
to within ±0.1 for temperatures ranging between 5 x 10^ K and 3x 10^ K and densities ranging 

between 10^ cm~^ and 10^° cm^'^ (HTimmer & Storey 1987). Our measurement is consistent 
with that of Soifcr ct al. (2004) who obtained a Spitzer 1R.S spectrum of the z = 3.91 quasar, 
APM08279+5255, and found Paa/Pa/3 = 1.05 ± 0.2 which is also too small for Case B 
recombination. 



5. Radio- loud versus Radio-quiet subsamples 

Twenty-two of the twenty-seven quasars in our sample are covered by the FIRST radio 
survey (Becker et al. 1995). Eight are detected, with 20 cm flux densities ranging from 
1.3 to 166 mJy; the remaining sources have upper limits 6*20 < 1-0 mJy The flve objects 
falling outside the FIRST survey area are covered by NVSS (Condon et al. 1998); one is 
detected with a flux density of 2.4 mJy (the remainder have S20 below this value). We 
compute the radio loudness parameter for each source starting with the formalism of Ivezic 
et al. (2002) which defines Rm — log(Fradio/-P'opticai) = 0.4(m — t), where m is any of the 
SDSS magnitudes, and t is the AB radio "magnitude" of the FIRST flux density, provided 
in the DRl catalog from which we selected our sources. We use the radio-loudness parameter 
without the logarithm (i.e., 10^™) and report its value in column (13) of Table 2. All have 
i? > 1.0. We do not include a /^-correction when computing Rr due to the unknown radio 
spectral indices for these quasars. Instead, we examined R^ for all flve SDSS fllters and 
found that in no object did R^n range over more than a factor of ~ 2. In our analysis we use 
Rr from the SDSS r band. 

Deflnitions for what constitutes a "radio-loud" quasar range from i? > 10 to i? > 100 
while radio-quiet quasars have R < 1.0. Depending on the deflnition, our sample contains 
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between two and five radio-loud objects. The rest of the sample has i? < 10 and are either 
considered radio- intermediate or radio-quiet. Although we have only a small sample of radio- 
loud objects, we can compare radio-loud and radio-quiet composites to look for any obvious 
differences. 

We have created two composites from the radio-loud (including all five objects with 
R > 10.0) and radio-quiet (including all other spectra) subsamples. The top panel of Figure 
9 plots the radio-quiet composite in blue with the radio-loud composite overplotted in red. 
We plot these in log(Luminosity x u) vs. log(z/) space which is standard for plotting spectral 
energy distributions. The radio-loud spectrum is scaled to the radio-quiet spectrum at 
1/ = 1.8 X 10^^ Hz, or 1.67/im, roughly where the signal-to-noise of our composite is highest. 

Superficially the spectra appear to show only small differences. This is likely a result of 
the small number of radio-loud quasars and our liberal definition of "radio-loud" {Rr > 10.0); 
a larger sample of radio-loud quasars, with Rr > 100 may reveal significant differences 
between the spectra. 

We compare our quasar composite spectrum made from the entire sample with the SED 
for quasars derived by Elvis et al. (1994) in the lower panel of Figure 9. For consistency 
the SEDs are scaled to our quasar composite at the same wavelength used in the top panel. 
The radio-loud SED is plotted with a red dotted line and the radio-quiet SED is displayed 
with a blue dashed line. Some of the differences observed between the radio-quiet and radio- 
loud subsamples in the top panel are present in the SEDs of the bottom panel, particularly 
below ~ 1.6 X 10^^ Hz (log(i/) < 14.2). The radio-loud quasars have a harder spectrum in 
the ~ 1.2 — 2 ^m band. Spectra obtained at longer wavelengths will determine whether this 
trend continues into the mid-infrared. In the SEDs of Elvis et al. (1994), when normalized in 
the same way as in Figure 9, the trend reverses around 10^'^ Hz, below which the radio-loud 
continuum gains strength out into the radio regime, as expected. Above the second break in 
the spectra, around 2.5 x 10^^ Hz (log(z/) > 14.4), the radio-quiet and radio-loud SEDs from 
Elvis et al. (1994) have harder continua than our composite continuum. The difference may 
be a consequence of the higher median luminosities of the quasars chosen for the Elvis et al. 
(1994) SED. Their radio-loud SED has 18 radio-loud quasars with a median My = -25.3 
and 29 radio-quiet quasars with a median My = —23.8. The median luminosity of our 
sources, in the equivalent SDSS g'-band is Mg — —23.4. 

We also fit the power-law plus blackbody model to our radio-loud and radio-quiet com- 
posites. We report these results in Table 5. The fit to the radio-quiet spectrum yields a 
softer spectral-index, a = —1.34, but a dust temperature consistent with that found for the 

total composite, Tdust = 1260 K. Since the radio-loud composite covers a smaller wavelength 
span and lacks coverage in the optical portion of the spectrum where the power-law domi- 
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nates, the fit to equation (5) does not yield physical results {a — —1.42, and Tdust = 5 K). 
Instead, we fit both composites to equation (5) fixing either the spectral-index or the dust 
temperature to the value derived for the total composite, a = —0.92 or T^ust = 1260 K. 
This approach allowed for fewer free parameters and generated resonable results, consistent 
with the total composite. When one parameter was held fixed, the radio-loud spectrum was 
fit by either a T^ust hotter by 100 K or a slightly softer spectral index compared with the 
radio-quiet spectrum. 



6. An Optical-to-Near-Infrcired Composite 

We used our composite-spectrum-building algorithm described in §3 with the optical 
spectra provided by SDSS and our IRTF spectra to build a composite quasar spectrum 
spanning the wavelength range 2700-37000A. We applied the same wavelength binning as in 
Equation (1) with Aq = 2700 A and n — 11250. The optical spectra were incorporated into 
the composite individually as separate spectra. They were first ordered by increasing redshift 
as for the infrared counterpart. Then, since each optical spectrum overlaps its infrared 
counterpart between roughly 8OOOA and 9500A, we used this region to scale and incorporate 
the first optical spectrum into the composite. After the full spectral range is filled, each 
additional spectrum (infrared or optical) is scaled to the previous partial composite using 
the median values of their mutually overlapping regions. 

Figure 10 shows the resulting spectrum plotted as log-wavelength versus fiux. We also 
plot a model host-galaxy composite spectrum (dashed line) to represent an estimated host 
galaxy contribution to our quasar composite. We generated the galaxy spectrum according 
to the following prescription: We assume that each quasar in our sample is hosted by a 
massive elliptical with a luminosity of 2 x . This is supported by observations of nearby 
quasar hosts (cf. Percival et al. 2001; Dunlop et al. 2003; Floyd et al. 2004). We use M^^ = 
—23.47 for early type galaxies (Kochanek et al. 2001) and scale the eUiptical galaxy template 
spectrum from Mannucci et al. (2001), which spans ultraviolet to near-infrared wavelengths, 
by its X-band flux to the flux expected from an 2 x L* galaxy at the redshift of each quasar. 
We then scale the final template by 20% to account for the fraction of galaxy light entering 
through our slit. At its maximum, the putative host galaxy contributes ~ 6% of the overall 
continuum around 0.9-1.0 fim (cf. §2). Table 7 lists the arithmetic and geometric mean 
spectra, available in their entirety in the electronic edition of the Journal. 
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6.1. K-corrections and Infrcired Colors 



A particularly useful application for this quasar composite is in calculating /^-corrections 
to passbands that lie in our spectral region. These include the three 2MASS bands, J (1.25 
/xm), H (1.65 /xm), and Ks (2.17 /xm), as well as the shortest Spitzer IRAC infrared band 
centered at 3.6//m. The long-wavelength baseline of our spectrum allows for X-corrections 
to be computed out to z ~ 3 — 6 for the 2MASS filters and z ~ 10 for the 3.6 //m band. 

Since our composite spectrum only extends to 3.5/xm and the IRAC S-G/im bandpass 
extends to ~ 4/xm, we must extrapolate the long wavelength end of our spectrum with the 
power-law fit for that region (derived in §4.1) to measure the rest- frame emission through 
the bandpass and relate redshifted emission to it. As there are no prominent emission lines 
expected in the 3.5-4 /im band, this continuum extrapolation should be adequate. 

We define the i^-correction, K{z), as the difference between the rest- frame absolute 
magnitude, M, of a source and its observed magnitude, m, at redshift z, through the same 
bandpass. Symbolically, 

m = M + 51ogio(^)+i^(^), (6) 

where Dl is the luminosity distance. The X-correction is then computed from the source 
flux per unit wavelength, /(A), and the response function for the bandpass, T(A) with the 
following equation: 



X(z) = 2.51ogio(l + ^) + 2.51ogio 



jT{X)f{X)dX 
jT{X)f(X/{l + z))dX 



(7) 



We use the response functions for the 2MASS bands and the Spitzer 3.6 fim IRAC band 
in equation (7) and compute the iiT-corrections for our composite quasar from 2; = — 10 in 
increments of Az = 0.1. The results are shown in Figure 11 and Table 8. 



7. Summciry and Future Work 



We have presented the first-ever near-infrared quasar composite spectrum spanning the 
wavelength range 0.57-3.5;um. Combined with optical spectra, this composite extends to 
0.3/xm. We have analyzed the continuum in this wavelength range to reveal three distinct 
power- laws with breaks around 1.08//m and 2.22/im. The optical part of the spectrum 
(blueward of 1.08 //m) is best fit by a spectral slope of ai, — —0.78 while the middle of the 
spectrum is best fit by — —1.81. For wavelengths longer than 2.22 //m, the spectrum is 
best fit by a^, — —1.03. We find that the order in which we combine our spectra to produce 
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the composite has the strongest effect on the optical spectral index, with the spectrum sorted 
from high to low rcdshift having a spectral index as hard as a^, = —0.37. The other spectral 
regions vary by Aa^ < 0.2. An alternative fit combining a single power-law and a blackbody 
provides an excellent description of the spectrum and offers strong evidence for the existence 
of hot dust (T(i ~ 1250 K) in these blue quasars. 

We have also studied the emission lines in this new spectral region, revealing the Paschen 
series lines as well as oxygen and helium and forbidden sulfur emission. Paschen line ratios 
rule out case B recombination in the broad-line region. The small number of spectra (13) 
obtained in the LXD band resulted in a modest signal-to-noise ratio beyond 2.1yum for our 
composite. Additional spectra in this band will allow measurements of the intensities and 
equivalent widths of the Brackett series lines expected in this wavelength region. 

We have compared radio-quiet and radio-loud subsample composites to the radio-quiet 
and radio-loud quasar SEDs produced by Elvis et al. (1994) . We find general agreement in 
the spectral shapes beyond ~ 1 fim. 

We have computed /^-corrections to the 2MASS J, H, and Kg bands as well as the 
Spitzer 3.6/im IRAC band for the purpose of predicting fluxes for high-redshift quasars 
observed in these bands. Since IRTF has an LXD mode which captures the 2.3-5.5/xm wave- 
length range, additional spectra in this band could extend our composite (for sources with 
the same redshift range) out to ~ 4.5//m. This will allow for a better measurement of the 
X-correction to the 3.6//m IRAC band as well as for the determination of the X-correction 
for the 4.5//m IRAC band when we extrapolate out to 5/im (the extent of the passband). 
Such a composite can be combined with future Spitzer IRS spectra of low-redshift luminous 
blue quasars planned by guaranteed time observers to extend the wavelength coverage, and 
our understanding of the intrinsic quasar spectrum, out to ~ 40/im. 
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A. Appendix 

Because near-infrared spectroscopy for the objects used to form the composite are not 
published in the literature, we make available the twenty-seven fully reduced quasar spectra 
to the astronomical community. Table Al displays a sample spectrum of SDSS J000943. 1—090839. 2. 
This spectrum, as well as the other 26 spectra used to create the composite, are available in 
their entirety in the electronic edition of the Journal. 
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Table 1. Observation Log 
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Table 2. Summary of Quasar Properties 
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17 


00 


16 


97 


16 


80 


16 


31 


16 


45 


14 


71 


13 


68 


12 


65 






—23 


15 


01 59 


10 


05 


+01 


05 


14 


5 





217 


17 


48 


17 


44 


17 


24 


16 


95 


17 


05 


15 


42 


14 


53 


13 


51 






—23 


16 


01 59 


50 


24 


+00 


23 


A A 

40 


8 





163 


15 


91 


15 


88 


15 


97 


15 


75 


15 


82 


1 A 

14 


07 


13 


09 


12 


01 


22.55 


15.2 


—23 


69 


02 17 


07 


87 


—08 


47 


43 


4 





292 


17 


78 


17 


41 


17 


14 


17 


13 


16 


85 


15 


77 


15 


08 


14 


15 






—23 


67 


02 42 


50 


85 


—07 


59 


14 


2 





378 


17 


23 


17 


06 


17 


03 


17 


09 


16 


71 


15 


88 


15 


12 


14 


16 






—24 


33 


03 12 


09 


20 


-08 


10 


13 


8 





265 


17 


72 


17 


51 


17 


32 


17 


24 


16 


87 


15 


98 


15 


16 


14 


17 






-23 


40 


03 22 


13 


89 


+00 


55 


13 


4 





185 


17 


08 


16 


86 


16 


61 


16 


14 


16 


34 


14 


81 


14 


01 


12 


90 






-23 


78 


03 40 


00 


50 


-05 


17 


47 


1 





217 


18 


12 


18 


00 


17 


76 


17 


76 


17 


35 


16 


23 


15 


39 


14 


22 






-23 


49 


07 36 


23 


12 


+39 


26 


17 


8 





118 


16 


26 


16 


39 


16 


25 


15 


49 


15 


77 


14 


18 


13 


33 


12 


27 


3.50 


3.0 


-23 


26 


12 55 


19 


69 


+01 


44 


12 


2 





343 


16 


02 


16 


01 


16 


05 


16 


25 


15 


52 


15 


23 


14 


47 


13 


43 


1.35 


1.0 


-24 


93 


13 07 


56 


57 


+01 


07 


09 


6 





276 


18 


10 


17 


58 


17 


15 


17 


06 


16 


50 


15 


44 


14 


67 


13 


57 






-23 


61 


15 06 


10 


50 


+02 


16 


49 


9 





135 


16 


47 


16 


39 


16 


32 


15 


85 


15 


98 


14 


68 


13 


72 


12 


79 






-23 


21 


16 05 


07 


92 


+48 


34 


22 








295 


16 


66 


16 


69 


16 


68 


16 


80 


16 


28 


15 


27 


14 


41 


13 


38 






-24 


01 


17 04 


41 


37 


+60 


44 


30 


5 





372 


15 


65 


15 


44 


15 


25 


15 


18 


14 


84 


14 


15 


13 


47 


12 


43 






-26 


19 


17 27 


11 


81 


+63 


22 


41 


8 





218 


17 


38 


17 


19 


17 


02 


16 


74 


16 


74 


15 


54 


14 


92 


14 


17 


1.36 


2.4 


-23 


37 


21 18 


43 


24 


-06 


36 


18 








328 


17 


65 


17 


67 


17 


59 


17 


69 


16 


80 


16 


38 


15 


65 


14 


45 


75.52 


225.1 


-23 


66 


21 26 


19 


65 


-06 


54 


08 


9 





418 


17 


16 


16 


92 


16 


85 


16 


76 


16 


62 


15 


94 


15 


27 


14 


39 


69.80 


105.6 


-25 


04 


23 22 


35 


66 


-09 


44 


38 


1 





372 


17 


11 


16 


99 


17 


05 


17 


05 


16 


55 


16 


12 


15 


23 


14 


19 






-24 


33 


23 49 


32 


77 


-00 


36 


45 


8 





279 


17 


38 


17 


25 


17 


12 


17 


20 


16 


84 


15 


97 


15 


28 


14 


02 






-23 


50 



Table 2 — Continued 



R.A. 


Dec 




u 


9 


r 


i 


z 


J 


H 


Ks 






M, 


(J2000.0) 


(J2000.0) 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag)*^ 


(mag)'' 


(mag)'' 


(mJy) 




(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


23 51 56.12 


-01 09 13.3 


0.174 


15.68 


15.76 


15.79 


15.37 


15.59 


14.33 


13.41 


12.18 


166.50 


94.6 


-24.23 



^Redshifts are from the SDSS quasar catalog (Schneider et al. 2003). 

''2MASS measurements are Vega-based, not AB, magnitudes. 
'^Fhix density from NVSS (Condon et al. 1998). 

'^Radio loudness parameter, defined by Rr = f radio / fr, where fr = 3631 x 10"° "^'' Jy. 
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Table 3. Composite Quasar Spectrum, Arithmetic and Geometric Means 



A 


<fx> 


< fx> Uncertainty 


< fx >gin 


(A) 


(arbitrary units) 


(arbitrary units) 


(arbitrary units) 


5900.4 


1.9824 


0.1895 


2.2860 


5901.7 


2.0507 


0.2737 


2.4589 


5903.1 


2.6955 


0.2954 


3.0849 


5904.5 


1.8391 


0.2292 


2.5981 


5905.9 


3.2544 


0.2552 


3.0030 


5907.2 


2.6429 


0.2391 


2.6205 


5908.6 


2.4293 


0.1401 


2.4708 


5910.0 


2.6892 


0.1674 


1.4671 


5911.4 


2.5463 


0.1675 


2.6436 


5912.7 


2.1509 


0.1073 


1.1448 



Note. — The complete version of this table is in the electronic edition 
of the Journal. The printed edition contains only a sample. 



Table 4. Fitted Power-Law Indices 



Spectrum 


UV - Optical 


a„ 


H-a 




mid-infrared 


ai, 


far-infrared 


a„ 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 



References. — (1) 



1 Lya - H/3 -0.46 ••• ••• ••• ••• ••• (4) This papci 

2 1350-4230A -0.44 6005 - 7180 A -2.45 • • • • • • • • • (6) This paper, near 

3 • • • • • • 6005 - 7180 A -1.03 ± 0.01 • • • • • • • • • (8) This paper, near 

4 ••• ••• 6005-7180 A -1.55 ±0.01 ••• ••• ••• 

5 ••• ••• 6005 - 7180 A -1.21 ±0.01 ••• ••• ••• 

6 ••• ••• 5700- 10850A -0.78 1.085-2.230 /xm -1.81 2.23-3.50 /xm -1.03 3.39 

7 ••• ••• 5700 -9570 A -0.37 0.957-2.093 /xm -1.73 2.09-3.50 /xm -0.96 3.39 

8 ••• ••• 5700 -9730 A -0.48 0.973-2.382 /xm -1.74 2.38-3.50 ;um -1.17 3.36 



I 

I 
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Table 5. Fitted Power-Law plus Blackbody 



Spectrum 




Tdust (K) 




(1) 


(2) 


(3) 


(4) 


1 


-0.92 


1260 


3.18 


2 


-0.94 


1260 


3.21 


3 


-0.93 


1250 


3.17 


4 


-1.34 


1260 


3.55 




-0.92^^ 


1230 


2.96 




-0.80 


1260* 


2.93 


5 


-0.92=* 


1330 


4.30 




-0.85 


1260* 


4.64 



*This parameter is held fixed in the fit. 

References. — (1) Near-infrared com- 
posite ordered by ziow — zuigh ; (2) Near- 
infrared composite ordered by Zhigh — 
ziow\ (3) Near- infrared composite or- 
dered by Zhigh snr - ziow snr, (4) Radio- 
quiet subsample; (5) Radio-loud subsam- 
ple. 
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Table 6. Emission Line Features 



Emission Feature 


Lab A 


Measured A 


Start'' 


End" 


Intensity 


EW 


FWHM 




(A) 


(A) 


(A) 


(A) 


[100 X F/F(Ha)] 


(A) 


(A) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Ha 


6563 


6569.5±0.3 


6463.6 


6674.5 


100.00±0.98 


180.4 


47.4 


01 


8446 


8457.5±0.3 


8228.7 


8684.9 


3.78±0.16 


10.7 


104.2 


[Q TTTl 

[b iiij 


yuoy 


9070. 8±0. 8 


nnc A 1 
9054.1 


9098.4 


0./5±0.04 


0.8 


inn 

ly.o 


Fell 


9202 


9214.0±1.4 


9115.4 


9312.5 


1.10±0.09 


3.5 


81.4 


Pae 


9545 


9534.4±0.5 


9443.4 


9627.5 


2.02±0.09 


7.0 


39.9 


PiiS 


10049 


10042.3±0.3 


9688.1 


10395.5 


5.72±0.13 


21.1 


161.2 


He I'' 


10830 


10830.0 


10635.2 


11025.6 


10.05±0.20 


36.0 


115.8 


Pa7^ 


10941 


10941.0 


10702.2 


11178.0 


1.95±0.71 


7.0 


109.4 


01 


11287 


11296.4±1.3 


11211.8 


11380.0 


0.89±0.04 


3.3 


78.8 


Pa/3 


12820 


12821.3±0.1 


12533.2 


13108.7 


4.82±0.06 


18.4 


128.3 


Paa 


18756 


18735.5±2.1 


18499.5 


18970.1 


3.10±0.01 


12.7 


196.1 


He I 


20580 


20506.8±6.4 


20035.8 


20975.2 


1.89±0.03 


8.6 


300.5 



''Wavelengths between which the line is integrated. The limits are chosen at 5% of the peak, above the 
continuum. 

''The line centers were fixed when fitting He I and Pa7 due to the strong blending of these lines. 



Table 7. Optical-to-near-infrared Composite Quasar Spectrum, Arithmetic and Geometric 

Means 



A 


< /a > 


< f\ > Uncertainty 


< fx >gm 


(A) 


(arbitrary units) 


(arbitrary units) 


(arbitrary units) 


2761.6 


12.2489 


0.3339 


8.9927 


2762.2 


11.9368 


0.3186 


10.7421 


2762.9 


12.0142 


0.3266 


10.7801 


2763.5 


12.5634 


0.3259 


11.1975 


2764.2 


12.4229 


0.3171 


11.0970 


2764.8 


11.8648 


0.3089 


10.7094 


2765.5 


11.9007 


0.3047 


10.6902 


2766.1 


11.9286 


0.3036 


10.6816 


2766.8 


11.9703 


0.3048 


10.7700 


2767.4 


12.3445 


0.3069 


11.0721 



Note. — The complete version of this table is in the electronic edition 
of the Journal. The printed edition contains only a sample. 
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Table 8. Quasar X-corrections 



Redshift 


K{J) 


K{H) 






0.1 


0.005 


0.085 


0.015 


-0.012 


0.2 


0.030 


0.111 


0.094 


0.001 


0.3 


0.085 


0.178 


0.202 


-0.005 


0.4 


0.085 


0.201 


0.279 


-0.004 


0.5 


0.071 


0.199 


0.320 


0.027 


0.6 


0.022 


0.210 


0.335 


0.050 


0.7 


-0.208 


0.294 


0.378 


0.065 


0.8 


-0.397 


0.287 


0.446 


0.095 


0.9 


-0.459 


0.275 


0.414 


0.136 


1.0 


-0.478 


0.260 


0.389 


0.190 


1.1 


-0.393 


0.231 


0.418 


0.253 


1.2 


-0.310 


0.196 


0.505 


0.299 


1.3 


-0.403 


-0.033 


0.516 


0.338 


1.4 


-0.466 


-0.262 


0.494 


0.372 


1.5 


-0.520 


-0.297 


0.485 


0.377 



Note. — The complete version of this table is in the 
electronic edition of the Journal. The printed edition 
contains only a sample. 
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Table Al. Spectrum of SDSSJ000943. 1-090839. 2 





A 


h 




fx Uncertainty 




(A) 


(arbitrary units) 


(arbitrary units) 


n 
u. 


,OU ( oou 


-1.07676c- 


-16 




0, 


,808083 


-3.07219c- 


-16 


1.19947C-16 


0. 


,808285 


-3.00444c- 


-16 


4.82051e-16 


0. 


,808488 


4.34456e- 


■16 


2.07032e-16 


0. 


,808690 


4.09549e- 


■16 


3.45445e-16 


0. 


,808893 


2.37477e- 


■16 


2.079976-16 


0. 


,809095 


3.07857e- 


■17 


1.896396-16 


0. 


,809298 


2.68657e- 


■17 


3.917886-16 


0. 


,809500 


-2.09553c- 


-17 


1.89665C-16 


0, 


,809702 


-6.93002c- 


-17 


9.67822C-17 


0. 


,809905 


-1.22391c- 


-16 


3.68280C-16 


0, 


,810107 


1.84263c- 


-16 


4.10956C-16 


0. 


,810310 


6.95212e- 


■16 


1.074316-15 


0. 


,810512 


-5.50372e- 


■17 


4.675146-16 


0. 


,810715 


-3.68670e- 


■17 


5.784696-17 



Note. — This spectrum, as well as the other 
26 spectra used to create the composite, are avail- 
able in their entirety in the electronic edition of 
the Journal. The printed edition contains only this 
sample. 
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1.0 1.5 2.0 2.5 3.0 3.5 

Wovelength (microns) 

Fig. 1. — Near- infrared quasar composite spectrum. The dotted lines show the positions of 

expected prominent emission Unes: Ha 6563, O I 8446, [S III] 9069, Fe II 9202, Pae 9545, 
Fa6 10049, He I 10830, Pa7 10941, O I 11287, Pa/3 12822, Paa 18756, He I 20520, Br7 21654, 
Br/3 26260 A. 
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0.20 




1.0 1.5 2.0 2.5 3.0 3.5 

Wavelength (microns) 



Fig. 2. — Top - Combined error spectrum. Bottom -Signal-to-noise per pixel of our com- 
posite quasar spectrum. 
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1.0 1.5 2.0 2.5 3.0 3.5 

Wavelength (microns) 

Fig. 3. — Near- infrared quasar composite spectrum, generated from the geometric mean of 
the individual spectra. The dotted hnes show the positions of expected prominent emission 
hues: Ha 6563, O I 8446, [S III] 9069, Fe II 9202, Pae 9545, P&S 10049, He 1 10830, Pa7 10941, 
O I 11287, Pa/3 12822,PaQ; 18756, He I 20520, Br7 21654, Br/3 26260 A. 




Wavelength (microns) 

Fig. 4. — Number of spectra contributing to the final composite as a function of wavelength. 
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3000 4000 5000 6000 7000 8000 

Wavelength (Angstroms) 



Fig. 5. — Upper - Optical composite spectrum generated from the 27 objects in the near- 
infrared sample (black). Overplotted are composites which we generated from the 87 SDSS- 
DRl quasars that obeyed our selection criteria (blue), and the SDSS quasar composite spec- 
trum from Vanden Berk et al. (2001) [red). Lower - Difference between the optical composite 
from this work (black line above) and the SDSS-DRl composite spectrum (red line above) 
(black). The template spectrum of a giant elliptical galaxy from Mannucci et al. (2001) is 
overplotted (with an arbitrary scaling, green) for comparison. The residual spectrum reveals 
both continuum and absorption features obvious in the galaxy spectrum. The absorption 
lines marked with a dotted line are Ca II K 3933, Ca II H 3968, CH G-band 4300, Fe I 4384, 
Mg I 5175, Na I 5892 A. The prominent TiO absorption between 7050-7275A is also marked 
by the shaded band. 
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Wavelength (microns) 

Fig. 6. — Composite quasar spectrum plotted in — X on logarithmic axes with power-law 
fits to three spectral regions overplotted in blue (0.58 — 0.97 /im), green (0.97 — 2.37 /im) 
and red (2.37 — 3.50 /im). The thick black lines at the bottom mark the continuum windows 
used in the power-law fit. 



0.6 0.7 0.8 0.9 1.0 2.0 

Wavelength (microns) 



3.0 



Fig. 7. — Power-law plus blackbody fit to our composite quasar spectrum plotted in F^, — A 
on logarithmic axes. The combined fit is shown in red, the power law component with a 
best-fit spectral index of a = —0.92 is shown in green and the blackbody with the best-fit 
dust temperature of = 1260 K is shown in blue. 
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0.63 0.64 0.65 0.66 0.67 0.68 
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1.24 1.26 1.28 1.30 1.32 





2.00 2.02 2.04 2.06 2.08 2.10 

Wavelength (/xm) 

Fig. 8. — Fitted emission-lines displayed over their fitted regions. The fitted profile is 
overplotted as a dashed line {red dotted line in electronic edition) and the continuum is 
plotted with a dash-dot line {blue dash-dot in electronic edition). In the He I-|-Pa7 panel, we 
overplot the individual line profiles with a dotted line. 
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Fig. 9. — Top - Comparison of composites made from radio-quiet (blue) and radio- loud 
{red) subsamples plotted in log(i/) vs. logi^uL^). Bottom - The geometric mean composite 
spectrum from the entire sample overplotted with the mean spectral energy distributions 
from Elvis et al. (1994) for radio-loud {red dotted line) and radio-quiet {blue dashed line) 
quasars. 
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Wavelength {fj-m) 



Fig. 10. — Optical-to-near-infrared quasar composite spectrum. Expected prominent emis- 
sion lines include {vertical /znes):[0 II] 3727, H5 4102, H7 4341, H/3 4862, [O III] 4959, 
[O III] 5007, Ha 6563, O I 8446, [S III] 9069, Fe II 9202, Pae 9545, PaS 10049, He I 10830, 
Pa7 10941, Pa/3 12822, Paa 18756, 617 21654, Br/3 26260 A. The dashed line shows the 
possible contribution from a 2 x L* eUiptical host galalxy. 
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Fig. 11. — i^T-corrections as a function of redshift for the three 2MASS bands and the Spitzer 
3.6//m IRAC band. 



